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ABSTRACT 


High bandwidth communications and optical/RF weapons systems are 
being developed that are limited by atmospheric absorption and accumulated 
phase distortions. The need and ability to mitigate these effects depends on the 
magnitudes of the atmospheric perturbations. The optical and RF index of 
refraction structure function parameters, Cn^, are examined simultaneously using 
field data. It is difficult to model the magnitudes of Cn, accurately and values 
produced by models are frequently off by an order of magnitude or more. In situ 
measurements of Cn^ are extremely sparse at RF wavelengths and limited to a 
few field studies. To refine models or conduct climatologically studies for Cn^ 
requires direct measurements to identify the underlying factors that produce the 
fluctuations and provide a clear understanding of the phenomena. This thesis 
utilized high speed measurements of the atmospheric humidity, temperature and 
wind speed collected on a tower 10 m above ground level at a coastal location. 
The humidity data were collected with a high speed infrared humidity sensor. A 
three axis sonic anemometer provided the wind fluctuation data and a fine wire 
temperature sensor as well as the sonic anemometer provided atmospheric 
temperature. All the data were sampled at 20 Hz. This study examined a subset 
of 251 days of data collected at the Marina, California airport to determine the 
relative variations of the optical and RF magnitudes of Cn^ and the underlying 
atmospheric phenomena that produced the results. 
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I. INTRODUCTION 


A. OBJECTIVES 

Understanding the refractive index structure function parameter (C^) and 

related parameters is essential when dealing with optical and radio propagation. 
A wide array of techniques for measurement of optical turbulence have been 
developed (Eaton 2005). High-speed point measurements of the temperature 
and humidity fluctuations provide a direct means to determine the optical (dry) 
and RF (wet) Cn^ simultaneously. 

The Naval Postgraduate School (NPS) Department of Meteorology, 
Monterey, CA, operates a sensor suite designed for flux measurements near the 
airport in Marina, CA. This suite includes a three-axis sonic anemometer, a fine- 
wire temperature probe and a humidity measuring system all taking 
measurements at a frequency of 20 Hz. The data collected from these sensors 
was used to measure and the related parameters directly. 

An important phenomena characterizing this and similar coastal overland 
locations is the land-breeze - sea-breeze wind reversal that occurs frequently, 
almost every day. Airflow data measured with the Marina site flux measurement 
suite from October 2005 to June 2006 were used to investigate how the optical 
and RF depends on the atmospheric wind flow and solar heating patterns. 

Furthermore, distinct sea breeze frontal passages were examined in detail to 
understand the abrupt changes in temperature and humidity that occur during 
these transitions. 
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B. BACKGROUND 


1. Index of Refraction 


The atmospheric index of refraction depends on pressure, temperature 
and humidity. The relative contribution of these variables changes with 
wavelength. For visible and near-infrared radiation, temperature fluctuations 
dominate the process. It is customary to use a reference wavelength value of 
0.5//m and the resulting equation for index of refraction is: 


^-1 = 79 x 10 ®^ 
T 


( 1 ) 


where nis the index of refraction, T is the temperature in degrees Kelvin, p is 
the pressure in mb (Beland 1996). Assuming that atmospheric turbulence occurs 
at speeds much less than the speed of sound, the change in the index of 
refraction n for an isobaric change in temperature is 


dn = 


r 

V 



y 


xlO^® 


( 2 ) 


At millimeter and centimeter wavelengths the effects of water vapor 
molecular rotations are substantial. Tatarski (1961) provided the following 
expression for the real part of the index of refraction that is valid for radio waves 
of centimeter wavelength: 


n-l = 10^'x 


79 


p + 


4800e^ 

T J 


(3) 


where n is the index of refraction, T is the temperature in degrees Kelvin, p is 
the pressure in mb and e is the water vapor pressure in mb. The corresponding 
change in index of refraction for changes in both temperature and the water 
vapor partial pressure is 


2 



dn = 


79 


p + 2'- 


379200 


dt + 


379200 


de 


X 10 


( 4 ) 


At RF wavelengths the humidity term often dominates the temperature term and 
the total change in n will depend on whether the changes in humidity and 
temperature add or subtract. 

2. Kolmogorov Turbulence Variance Statistics 

For many geophysics processes, the conventional variance is not a robust 
statistic because the mean value is not stationary. Kolmogorov introduced a 
more advanced statistic based on structure functions. By assuming homogeneity 
and isotropy at least in a local volume, and if the random processes vary slowly, 
structure functions represent the mean square difference in the fluctuations of 
f(ri,r 2 ) over a distance and r 2 between the two measurements (Tatarski 1961). 

D/(ri,r2) = ([/(ri)-/(r2)f) (5) 

Structure functions avoid the nonstationary mean issue and provide robust 
statistical descriptions of turbulent phenomena. Although the wind velocity is the 
dominant physical process involved in turbulence, the wind can carry passive 
additives like temperature and humidity with the turbulent flow. The passive 
additives pick up the same statistical characteristics from the turbulent flow and 
share similar statistical characteristics (Tatarski 1961). 

According to Kolmogorov turbulence theory, turbulent eddies range in size 
from macro-scale to micro-scale forming a continuum of decreasing eddy sizes. 
Energy from convection and wind shear is added to the system at the outer scale 
Lq (10’s - 100’s of meters) before it cascades to a smaller scale 4 (1cm or less) 
where viscosity converts the energy to heat. Fig. 1 (Max 2006). 
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Figure 1. Illustration of Kolmogorov turbulent eddies cascading down in the 
inertial dissipation range between the outer scale and the inner scale. 

(Max 2006) 


With dimensional arguments and assuming a homogeneous 
incompressible isotropic, medium, Kolmogorov showed that the longitudinal 
structure function of the velocity is: 

DM) = Cy'\ l,<r<y ( 6 ) 

cl is the structure function parameter that applies to the property specified by y, 

and r = \r^-rl\ is the spatial separation between the two measurements. The 

proportionality of the structure function in the inertial range {i^<r<LQ) also 

applies to the passive additive structure functions such as those for temperature, 
for humidity, and for refractive index (Tatarski 1961) 

3. Optical Turbulence 

Optical turbulence is the result of temporal and spatial fluctuations of the 
index of refraction resulting from atmospheric turbulence. In a basic differential 
form cl is a function of three elements and by squaring Eq. (4) is 


4 



( 7 ) 


^ dn^, 

2 

( dn^ 

2 

( dn^ 

'dn'" 


Cj + 


+2 



Wj 

1 

Idqj 

^ 1 

W) 

ydq, 


where represents the temperature structure parameter, C\ represents the 
structure parameter for humidity, and cross covariance between temperature 

and humidity. Eq. (7) is often written in terms of the potential temperature or 
virtual potential temperature particularly when considering variations over altitude 
such as with radar data (Stankov 2003). At optical frequencies, humidity 
generally has a negligible contribution except over saturated soil and maritime 
surfaces where its contribution to total C] can reach 20%, mainly through the 

term (Beland 1996). 

For Kolmogorov turbulence the refractive turbulence structure function 
parameter is 

= ( 8 ) 

which is the mean-square statistical average of the difference in the indices of 
refraction, n^andn^, between two points (Parker 2002). The angle brackets 

represent the ensemble average and r is the distance between the two points. 
Differences in the index of refraction are caused by localized changes in 
temperature and humidity resulting from turbulent mixing of the atmosphere. A 
typical value of r used to determine cl is one meter. 

cl is computed in the same manner yielding: 

Cf = ((r,-r,)")/r“ (9) 

By taking the partial derivative of the Eq. (1) with respect to the temperature and 
assuming the turbulent eddies are isobaric cl is related to the temperature 

structure function (C^) with the following equation: 
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(10) 


Cl={19x\Q-^Pnycl , 

where P is the atmospheric pressure in millibars and T is the temperature in 
Kelvin (Tatarski 1961). 

With high resolution temperature and humidity data measured 
simultaneously Eqs. (1) and (3) provide a means to compute Cn^ directly from 

c,;(11) 

This technique allows a direct comparison of the optical (dry) and RF (moist) Cn^ 
variations throughout a day as the underlying turbulent, atmospheric processes 
evolve by choosing the appropriate dry and moist expressions for n, Eqs (1) and 
(3) respectively. Part II of this thesis will introduce a more recent and refined 
version of the RF refractivity equation that was used for the data reduction. 


C. MEASUREMENT CHALLENGES 

1. Taylor’s Hypothesis 

Turbulence affects optical propagation through temporal and spatial 
distortion of the wave-front. Taylor’s hypothesis is an approximation that allows 
for estimating spatial turbulence statistics using measurements from air moving 
past a single point. For low intensity turbulence in a uniform mean flow the 
turbulence pattern is assumed to be “frozen”. This hypothesis is valid under most 
conditions, but becomes questionable when free convection velocity fluctuations 
are of the same order or larger than the mean velocity, or when the frequency of 
the turbulent eddies is not significantly higher than the magnitudes of the mean 
wind shear (Arya 2001). 
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2 . 


Nyquist Frequency 


The data sample rate fs determines highest frequencies that are resolved 
in a continuous time series. The Nyquist frequency fn is the highest frequency 
resolved in the sampled data set. To prevent aliasing a wave must be sample at 


least twice per period so the Nyquist frequency isf„ = y. 


3. Sampling Requirements 

Combining Taylor’s hypothesis guidance and the Nyquist frequency 
constraints poses a significant challenge when trying to measure turbulent eddies 
that range in size from centimeters up to meters with measurement at a single 
point. As the mean wind speed increases the required sampling frequency 
increases. Likewise, as the eddy sizes decrease the sampling frequency must 
increase. 

Satisfactorily resolving features in the centimeter range with typical mean 
flows of several meters per second requires instrument response times in 
milliseconds. Typically fine wire temperature probes, high frequency sonic 
systems, radars and lasers are used to measure such features at the necessary 
frequency. 

Table 1 shows the required sample frequency based on eddy wavelength 
and mean wind speed. The data used for this thesis was sampled at 20 Hz so 
frequencies in the highlighted areas represent what can be resolved from the 
data. 
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Required Sample Frequency (FIz) 


Mean flow 
U(m/s) 

Edc^ wavelength 

0.01 

0.05 

0.1 

0.25 

0.5 

1 

2 

0.5 

100 

20 

10 

4 

2 

1 

0.5 

1 

200 

40 

20 

8 

4 

2 

1 

1.5 

300 

60 

30 

12 

6 

3 

1.5 

2 

400 

80 

40 

16 

8 

4 

2 

2.5 

500 

100 

50 

20 

10 

5 

2.5 

3 

600 

120 

60 

24 

12 

6 

3 

3.5 

700 

140 

70 

28 

14 

7 

3.5 

4 

800 

160 

80 

32 

16 

8 

4 

4.5 

900 

180 

90 

36 

18 

9 

4.5 

5 

1000 

200 

100 

40 

20 

10 

5 

5.5 

1100 

220 

110 

44 

22 

11 

5.5 

6 

1200 

240 

120 

48 

24 

12 

6 

6.5 

1300 

260 

130 

52 

26 

13 

6.5 

7 

1400 

280 

140 

56 

28 

14 

7 

7.5 

1500 

300 

150 

60 

30 

15 

7.5 

8 

1600 

320 

160 

64 

32 

16 

8 

8.5 

1700 

340 

170 

68 

34 

17 

8.5 

9 

1800 

360 

180 

72 

36 

18 

9 

9.5 

1900 

380 

190 

76 

38 

19 

9.5 

10 

2000 

400 

200 

80 

40 

20 
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Table 1. Required sample rates as a function of mean wind and feature size. 

D. MARINA FIELD SITE SENSORS 
1. Location 

The NPS Marina Field Site sensor site is located 1.5 km north-north east 
of the runway at the Marina airport (formerly Fritzsche Army airfield). The sensor 
is approximately 4.7 km east of Monterey Bay and 40 meters west of a small bluff 
overlooking the Salinas valley, Fig. 2. The site is approximately 51 m above 
mean sea level and includes a 10 m tower for the flux measurements sensor 
suite. 
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Figure 2. Location of NFS sensor site in Marina CA. 

Wind at the site is primarily driven by two factors: Land and sea breezes, 
which have due east or westerly winds (Fort Ord TRFN 1977) and the orientation 
of the Salinas Valley, which produces north-east and south-west winds. The 
three sensors that collected the data for the calculations were mounted about 

9.5 meters above the surface facing a direction of 210°. This allows for the for 
minimal flow disturbance from the tower and support apparatus during the most 
common wind flows. Fig. 3. 
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Figure 3. 3-axis sonic anemometer with a fine wire temperature probe extending 
from the left side of the mount and the open path gas analyzer mounted to 
the right, (picture by Richard Lind of NPS) 

2. Sensors 

a. Ultra-Sonic Anemometer 

Sensors for examining Cn^ have to sample the influencing airflow 
properties fast enough to describe the relevant turbulent intensity and correlation. 
The sensors installed at the Marina Field site included a three-axis sonic 
anemometer, a CSAT3 manufactured by Campbell Scientific Inc. It uses three 
pairs of non-orthogonally oriented ultrasonic transducers. Horizontal wind velocity 
and direction along with vertical velocities are derived from the measured time of 
flight for the ultrasonic signals (manufacture’s instruction manual). 

The anemometer is capable of sampling at frequencies from 1 to 
60 Hz. During the measurements for the period studied, the anemometer 
sampled at 60 Hz and the data were block-averaged to achieve 20 Hz output 
data. Because the speed of sound is primarily dependant on temperature, the 
CSAT3 measurements also provided estimates of air temperature from the path- 
average speed of sound along the three sensor paths. This derived temperature 
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is density-based and referred to as the sonic virtual temperature (TJ, dependent 
of water vapor as well as temperature. The CSAT3 computes the sonic 
temperature, T^, in degrees Celsius, without this correction, with the following 
formula: 

T -273.15 (12) 

raK 

where c is the speed of sound, is the ratio of specific heat of dry air at 
constant pressure to that at constant volume (1.4), and is the gas constant 
for dry air (287.04 JK ' kg '). 

b. Fine Wire Thermocoupie 

The sonic anemometer has limitations with regard to evaluating 
optical turbulence from the influence of humidity on temperature estimations and 
from the measurement being a path average. This introduces a temperature 
uncertainty and a low-pass filter effect on the temperature turbulence statistics. 
A high response fine wire thermocouple, a Campbell Scientific FW1 collocated 
with the CSAT3 at the Marina Field site, also provided temperature fluctuations. 
The FW1 uses a 0.0254 mm wire. The fine wire’s small diameter allows 
millisecond response times and measurement at 20 Hz (Roper 1992). The small 
diameter also minimizes solar loading to the point where a shelter is not needed 
for the sensor (manufacture’s brochure), although the fine wire is susceptible to 
breakage. 


c. Open Path Gas Anaiyzer 

Because humidity fluctuations influence Cn^ and are a critical 
component of the RF atmospheric index of refraction, a high-speed sensor for 
humidity was important. The LI-7500 manufactured by Licor is an open path gas 
analyzer designed to measure Hp and CO^ in the air. It transmits an infrared 
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beam and measures attenuation as 4 wavelengths; attenuation at non-absorbing 
wavelengths (3.95//m and 2.40//m) is measured as a reference. Then 
absorption centered at 4.26 /um 2.59 /um is measured to determine CO^ and 

water vapor respectively (manufacture’s description). Measurements can be 
made a 5, 10, or 20 Hz, data collected for this project was 20 Hz. 

The LI-7500 also has a barometer housed in a junction box 
attached to the tower. This provides the pressure data used in this project. This 
sensor does not sample at high frequency. The scale of turbulence studied, the 
inertial subrange, is isobaric in nature so the high frequency pressure data is not 
required. 
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II. DATA AND ANALYSIS 


A. INTRODUCTION 

Understanding phenomena being examined and the capabilities and limits 
of sensors providing data requires a strategy for analyzing the data. Several 
hypotheses and assumptions are applied in the strategies. For example, using 
Taylor’s Hypothesis allows to be extracted by direct application of the 

structure function without using the more complex formulae needing power 
density spectra obtained from time series of varying signals. 

B. TURBULENCE EXTRACTION STRATEGY 
1. Application Taylor’s Hypothesis 

The first requirement for analyses of the fixed point measurement is an 
estimate of the size of the turbulent eddies causing the fixed point measured 
fluctuations. For the purposes of this study size a separation distance of 1.0 m 
was used for most of the data reduction. With a 20 Hz sample rate the number 
of sample points that make up a one meter segment is 

SP = = (13) 

U 

where SP is the number of sample points required, / is the sampling frequency, 
and U is the mean horizontal wind velocity. 

To determine U, a second order Butterworth filter was used with a 
cut of frequency of 0.001 of the Nyquist frequency. The low bypass Butterworth 
was used instead of simple averaging to limit the impact of extreme outliers and 
erroneous data. The U values obtained from this method are approximately 
equal to a one minute running average. 
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2 . 


Computing Refractivity (N) 


a. Radio Frequency 


For this work, Rueger provides an expression for the refractivity N 
that is more accurate than Eq. (3) (Rueger 2002): 


N = (n-l)xl0" 


77.689xP, 71.2925xe 375463xe 

-^ +^- 

rp rp 


(14) 


where is the partial pressure of dry air in millibars, T is the temperature in 

degrees Kelvin, and e is the partial air pressure of water vapor. This expression 
incorporates the current atmospheric CO 2 increases. 


b. Optical Wavelengths 


For optical wavelengths we simply use the dry component of 
Rueger’s radio frequency formula. 


N= 


77.689 X P^ 


(15) 


where P^ is the partial pressure of dry air in millibars, T is the temperature in 
degrees Kelvin. 


3. Computing 


a. Data and Point Selection 

derived from Eq. (8) requires a length within the inertial 

subrange, 1 m for this research, and values of refractivity from the two points at 
ends of that length. The main challenge is determining how many sample points 
make up that 1 m length. A secondary challenge is determining which 
temperature values to use. After these are done satisfactorily it is simply a 
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process of then applying the resulting data to the structure function parameter 
calculation, Eq. (8). 


b. Sonic Versus Fine Wire Temperatures 

A decision was made to use the sonic virtual temperatures values, 
obtained from the CSAT3, based on the discovered effects of solar radiation on 
the fine wire thermocouple. During hours of darkness and light winds the fine 
wire and sonic temperatures are nearly equal as would be expected based on 
theory the pressure temperature and humidity levels on the days examined. 
However after sunrise a noticeable difference arises between the two values. 
Fig. 4. In addition to the time of day associated with the difference, the 
conclusion that the difference is solar induced is further supported by the fact the 
difference nearly disappears after the sea breeze front passes and wind speeds 
increase. The stronger wind speed ventilates the fine wire sensors sufficiently 
and minimized the solar heating effect. 



Figure 4. Fine wire T (green) and sonic T (blue) for 1400-2300 UTC 

17 November 2005 
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The structure function parameter that defines is based on the 

spatial differences of passive variables. The spatial offset between the fine wire 
thermocouple and the measurement axes of the sonic anemometer was one 
factor in the decision not to use the fine wire data for determination. In 

addition, probe breakage reduced the available data to process. However the 
fine wire data provided a good way to verify the performance of the CASTS. The 
sonic temperature (CASTS) is derived from same sonic measurement used for 
wind velocity used for application of Taylor’s hypothesis, so problems with spatial 
and temporal offset do not arise. The issues with the sonic temperature, as 
discussed, are the influence of humidity fluctuations on the basic measurement 
and low-pass filter effect due to the path average. It should also be noted that the 
sonic derived temperature is altered by changes in the absolute humidity. 
However this difference is negligible in this situation. Fig. 5 illustrates the 
variations in temperature and humidity differences. 



0 5 10 15 20 25 30 

Time (minutes) 



Figure 5. Top: The percent change in the humidity (green) compared to the 

change in sonic T (blue) 

Bottom: The percent change in the fine wire T (red) compared to the 
change in sonic T (blue) for 2100-2130 UTC 23 February 2006 
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c. 


VERIFICATION OF Cl DETERMINATION STRATEGY 


1. Values 

The initial verification of the chosen extraction strategy was to check if 

the resulting values are within the range of values measured by other 

methods (10 to 10 " The results are consistent with expected values, as 
shown in the Fig. 6 example set (Walters 1991). 



Figure 6. Wet (green) and dry (blue) Cn^ values derived directly from the 
structure function, 2100-2115 UTC on 23 February 2006 

2. Power Spectral Density 

The values of are determined by passive variables (temperature and 

humidity) being mixed by turbulent flow. In the inertial sub range this turbulence 
follows the Kolmogorov theory. Therefore the log plot of power spectral density 
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for passive variables, or the quantities such as refractivity (N) determined by 
these variables, should show a dependence in the frequency domain. This 
was seen consistently as the example in Fig. 7 illustrates. 



Frequency (Hz) 

Figure 7. Power spectral density of Cn^ generate from 4096 data points (~3.5 
min of data) beginning at 2147 UTC 17 November 05 


During the time period corresponding to Fig. 7 the mean wind velocity was 
approximately 4 m s'^ for which, with the structure function distance set a 1 m, 
the resulting Nyquist frequency is 8 Hz. The departure from a slope at high 
frequencies, above 2 Hz, should not be interpreted as a departure from 
Kolmogorov theory. It is likely the result of noise in the signal as the sensor 
limitation (Nyquist frequency) is approached. The height of the data collection 
from the ground (9.8m) also imposes an outer scale at low frequencies. 
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D. TEMPERATURE AND HUMIDITY CORRELATION 

1. Radio Frequency and Optical Differences 

There were times when the RF was larger than the optical case. 

These events occurred with a negative correlation between temperature and 
humidity. Figures 8 and 9 provide a clearer illustration of these phenomena 
around the 45 minute point. When the optical is higher than the RF case the 
temperature and humidity correlation is positive. 



Figure 8. RF (green) and optical (red) from 2140-2155 UTC on 6 

December 2005 
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Figure 9. Temperature and humidity correlation coefficient from 2140-2155 UTC 

on 6 December 2005 

To further investigate the phenomena of the RF versus optical 

difference, the refractivity N 2 -N 1 in the moist and dry components of N 1 m apart 
were plotted. The last two terms of Eq. (14) are the moist component of N. The 
first term of Eq. (14) is the dry term (and also the only term in Eq. (15) for optical 
N). The difference of these components over one meter is plotted in Fig 10. 

0.15 r 


0.1 



40 45 50 55 

Time (minutes) 


Figure 10. Refractivity difference N 2 -Niof the separate dry N (blue) and moist N 
contributions to the refractivity (magenta) from 2140-2155 UTC on 6 

December 2005 
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2. Phase and Magnitude 


Fig. 8 and 9 show that when the optical exceeds the RF the moist 

and dry components of N are out of phase. This can only occur when 
temperature and humidity are positively correlated. The magnitude of the change 
humidity term in Eq. (14) is great enough to dominate the effects of the change in 
temperature and cause the moist and dry components of N to cancel. 

The temperature and humidity values between the 44 and 46 minute 
points are primarily negatively correlated (Figure 9), or the humidity fluctuations 
overwhelm the temperature changes. The influence of the humidity and 
temperature terms in Eq. (14) explains the large increase in the RF during 
these times. 
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III. CASE STUDIES 


A. DAY SELECTION 

The available multi-sensor data set spanned 251 days. Numerous days 
had missing or unreliable data for different sensors required to determine . 

These missing or unreliable data were caused by broken fine wires or 
contaminated lenses on the LI-7500. Within the set of days with all sensors 
operating, there were six that had clear skies and distinct sea breeze fronts. 
From these, two days (17 November 2005, and 23 February 2006) with what 
were interpreted to be the most well-defined see breeze fronts were selected for 
detailed analyses and interpretation of Cn^ variations and impacting factors. 

Neither of the previously mentioned days had a complete 24 hours of 
data. To view a diurnal cycle, a day with 24 hours of good data was analyzed (26 
November 05). 

B. 17 NOVEMBER 2005 

1. General Conditions 

The meteorological conditions on 17 November 2005 were nearly perfect 
for examining a sea breeze front Fig. 11. There were weak synoptic forcing 
conditions with offshore winds less than 4 m s \ Clear skies as indicated by the 
shortwave irradiance curve permitted strong solar heating as indicated on the 
temperature plot. At approximately 2145 UTC (1345 PST) the passage of the 
sea breeze front was clearly marked by a rapid change in wind direction from 
offshore (-080°) to onshore at the field site location (-270°) and by wind speed 
increasing (from less than 2 m s'^ to above 4 m s'^), and with increasing dew 
points and decreasing temperatures. Increases in the wind speed forced 
corresponding changes in Cn^ by causing increases in mixing of airflow 
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temperature and humidity resulting in lower vertical gradients in the overlying 
airflow. 
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Meteorology; Fort Ord (Profiler Site), CA (elev. 51 

Data courtesy Department of Meteorology, Naval Postgraduate School 
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<-Time Increases io Left 

Data from 17-N0V-2005 00:00 through 17-NOV-2005 23:58 UTC 
_ Image generated by Department of Meteorology, Naval Postgraduate School _ 

Figure 11. Weather summary for 17 November 2005 


2 . 


Sea Breeze Shift 


Examination of the 30-minute block containing the sea breeze the data 
from the CSAT-3 and the LI-7500 also shows the very distinct sea breeze as 
seen in Fig. 12. Over a period of approximately 30 seconds there is a 
temperature decrease (~ 4°K), an absolute humidity increase (~ 2 g m'^), and a 
doubling of the wind velocity from ~2 m s'^ to ~4 m s \ 
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Figure 12. Wind humidity and temperate from 2130-2200 UTC 


The time series of from the same period reveals two interesting 
features, Fig. 13. First, the optical values are actually lower than the 
RF values prior to the sea breeze. Second there is a massive spike in 
the RF values (~1.5 orders of magnitude) at the sea breeze front. 
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Figure 13. RF (green) and optical (red) from 2130-2200 UTC 
Eqs. (5), (11) and (12) for show that the RF includes the optical 

(solely temperature dependant) plus terms that take into account moisture with 
the moisture term in the numerator and the temperature values in the 
denominator. To achieve RF values lower than the optical C^, temperature 
and humidity values must be positively correlated and the changes in humidity 
small enough between sample points that the resulting change in RFC^ is less 

than the change in optical C^. If the temperature and humidity changes are 
negatively correlated the cross correlation term (term 3 in Eq. (7)) has a large 
value and produces RF refractivity fluctuations that are greater than seen for 
the optical case 

Fig. 14 is a plot of the temperature humidity correlation during 30 minutes, 
15 minutes before, during and 15 minutes after the sea breeze. Based on 
Kolmogorov theory the correlation coefficient will always be ± 1 (Hill 1989). 
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Because the correlation coefficients include an average of 30 seconds of data 
and the LI-7500 is ~10cm offset from the CSAT-3 paths, the correlation 
coefficients of exactly one are not shown. The overall change in the sign of the 
correlation coefficients is evident in the data. 



Figure 14. Temperature and humidity correlation coefficient from 2130-2200 UTC 

3. Sea Breeze Boundary 

The high frequency data allowed for a detailed examination of the frontal 
passage which lasted ~30 seconds. Within this 30-second period the power 
spectral density of N continued to follow the dependence expected by 
Kolmogorov theory, Fig. 15. The temperature and humidity correction coefficients 
averaged for every second showed a distinct shift at 45.4 minutes. Fig. 16. 
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Figure 15. Power spectral density of N as sea breeze boundary passed, 30 

seconds in the 2145Z minute 



Figure 16. Temperature and humidity correlation plotted every second during 

boundary passage. 
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c. 


23 FEBRUARY 2006 


1. General Conditions 


This day also had a sharp sea breeze front as seen in Fig. 17. There were 
weak synoptic forcing conditions with offshore winds less than 4 m s \ Clear 
skies as indicated by the shortwave irradiance curve permitted strong solar 
heating as indicated on the temperature plot. At approximately 2115 UTC (1315 
PST) the passage of the sea breeze front was clearly marked by a rapid change 
in wind direction from offshore (~080°) to onshore (~270°) at the field site location 
and by velocity increasing (from ~2 m s'^ to above 6 m s"^), and with increasing 
dew points and decreasing temperatures. As seen for 17 November 2005, the 
significant events that produced changes in were the wind speed increase 

that increased the mixing, and the airflow temperature and humidity changes 
causing different vertical gradients in the overlying airflow. 


Surface Meteorology: Fort Ord (Profiler Site), CA (elev. 51 m) 

Data courtesy Department of Meteorology, Naval Postgraduate School 
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Hour: 00232221 20 19 18 17 1 6 15 14 13 12 1 1 1009 0807 060504030201 00 

<-Time Increases to Left 

Data from 23-FEB-200B 00:00 through 23-FEB-2006 23:58 UTC 
Image generated by Department of Meteorology, Naval Postgroduate School 

Figure 17. Weather summary for 23 February 2006 


29 




























2 . 


Sea Breeze Shift 


The temperature shift is not as distinct and in the pervious case dropping 
~4°K in ~7 minutes. The winds jump from a mean value near 2 m s'^ to a main 
value near 7 m s'^ in the same 7-minute period. The absolute humidity values to 
make the sharp jump noted the previous case going from 9 g m-^ to 12 g m-^ in 
~30 seconds, Fig. 18. 
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Figure 18. Wind humidity and temp 2100-2130 UTC 
During this event the RF exceeded the optical before and after the 

sea breeze front passed, Fig. 19, with a noticeable spike (~1 order of magnitude) 
during the passage of the front. This also corresponds with a temporary reversal 
of the temperature humidity correlation coefficients during the passage. Fig. 20. 
Fig. 20 also shows that the correlation between temperature and humidity was 
mostly positive. Since Fig. 19 shows that the RF exceeded the optical, this 

implies that the humidity fluctuations had to exceed the temperature fluctuations 
throughout the entire period. 
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Figure 19. RF (green) and optical (blue) from 2100-2130 UTC 



Figure 20. Temperature and humidity correlation coefficient from 2100-2130 UTC 
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D. 


26 NOVEMBER 2005 


1. General Conditions 


This day had a stronger synoptic patter the resulted in westerly flow nearly 
all day. In the three hours prior to sunrise and the three hours after sunrise and 
easterly land set up for brief periods but was not able to maintain. The longwave 
irradiance plot indices some brief periods with clouds during the night, but skies 
were predominately clear during the period, Fig. 21. 


Surface Meteorology: Fort Ord (Profiler Site), CA (elev. 51 m) 

Data courtesy Department of Meteorology. Naval Postgraduate School 
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Data from 26-N0V-2005 00:00 through 26-NOV-2005 23:58 UTC 
Image generated by Department of Meteorology, Naval Postgraduate School 


Figure 21. Weather summary for 26 November 2005 
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2. 


Cn^ and Temperature-humidity Correiation 


During the overnight hours, the RF exceeded the optical by nearly 

one order of magnitude, Fig. 22, while the temperature and humidity correlation 
coefficients remained primarily negative. Fig. 23. In the three hours prior to 
sunrise, when the land breeze intermittently set up, correlations varied between 
positive and negative as the optical and RFC^ had similar magnitudes. After 
sunrise (~1600-1800 UTC) the winds remained light (< 2 m s^) and convection 
dominated the mixing processes. During that period RF increased by nearly 

2 orders of magnitude when the temperature and humidity were negatively 
correlated. 



Figure 22. 24-hour plot of , red: optical, green: RF 
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Figure 23. 24-hour plot of temperature and humidity correlations (90 sec 

averages) 
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IV. CONCLUSIONS AND RECCOMENDATIONS 


A. DATA AND EQUIPMENT 

This thesis studied the interaction of the atmospheric temperature and 
moisture fluctuations and their effects on the optical (dry) and RF (moist) 
refractive index structure function parameter, Cn^, using data obtained with a 
turbulence instrumentation suite at the NPS Marina Field Site, Marina Airport, 
Marina, CA. The instrumentation consisted of a sonic anemometer, fine wire 
thermocouple and open path gas water vapor analyzer sampled at 20 Hz., The 
equipment and sampling approach were suitable for determining values of 

directly from the structure function for both the dry and moist indices of refraction. 
The onset of the sea breeze produced a strong increase of one to two orders of 
magnitude in the RF refractive fluctuations, when the cross correlation between 
temperature and humidity changed signs. 

B. RECOMMENDATIONS FOR FUTURE STUDY 

This study, although limited in scope, indicates that further data 
processing could reveal the relationships between the underlying atmospheric 
conditions that produce the optical and RF refractivity fluctuations Since the 
correlation of temperature and humidity values play a significant role in 
enhancing the RF , a more extensive study could allow one to predict the sign 

of the correlation coefficient in advance. Identifying a larger number of days with 
data sets uninterrupted by dropouts or equipment failure, would allow the 
processing of more days. This could reveal additional diurnal and seasonal 
patterns in the data as well as the creation of meaningful averages. Data 
collected at other locations such as the East coast of the US would be beneficial 
since the inversion conditions and moisture concentrations are appreciably 
different than for the Pacific Coast. 
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